Effects of particle-size ratio on jamming of binary mixtures 
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We perform a systematic numerical study of the effects of the particle-size ratio R > 1 on the 
properties of jammed binary mixtures. We find that changing R does not qualitatively affect the 
critical scaling of the pressure and coordination number with the compression near the jamming 
transition, but the critical volume fraction at the jamming transition varies with R. Moreover, the 
static structure factor (density correlation) S(k) strongly depends on 7? and shows distinct long 
wave-length behaviors between large and small particles. Thus the previously reported behavior of 
S(k) ~ k in the long wave-length limit is only a special case in the R — > 1 limit, and cannot be 
simply generalized to jammed systems with R > 1. 
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I. INTRODUCTION 

At zero temperature and shear stress, a packing of fric- 
tionless spheres undergoes the jamming transition when 
its volume fraction <f> is varied across a critical value <j) c 
This transition is denoted as Point J in the jam- 
ming phase diagram 2, 3]. Although c/) c may vary with 
protocols [1, @ , it has been demonstrated that the max- 
imally random jammed packings of hard spheres with 
single size (mono-disperse) usually exist at a well-defined 
volume fraction (f> c 0.64 which is also the most 
probable volume fraction for hard spheres to jam at [6j. 
When the spheres are deformable and interact via pair- 
wise short-range repulsion, a marginally jammed solid 
is formed when the volume fraction is slightly greater 
than <j) c . Marginally jammed solids have multiple critical 
scalings with (f> — <p c oi some mechanical and geometri- 
cal quantities such as the pressure, bulk modulus, shear 
modulus, and coordination number (average number of 
interacting neighbors per particle) [lj, implying the crit- 
icality of Point J. Point J is isostatic [l|, i.e. the coor- 
dination number z — z c — 2d, where d is the dimension 
of space. The isostaticity is the minimum requirement 
of particle constraints to maintain mechanical stability 
and plays a key role in the control ofproperties of typi- 
cal jammed systems such as glasses 0, |H,[!|. Therefore, 
understanding Point J is a good starting point to study 
some long-standing problems, e.g. glass transition [10| 
and the unusual thermal properties of glasses [ll| . 

The jamming transition at Point J is puzzling com- 
pared to conventional critical phenomena. For instance, 
during the jamming transition, there is no apparent 
structural change. Consequently, a diverging static corre- 
lation length, e.g. the density correlation length present 
in the long wave-length limit, that can be observed in 
typical phase transitions is lacking. However, jammed 
packings of mono-disperse spheres still exhibit nonan- 
alytic density fluctuations different from simple liquids 
[12, HH: the static structure factor S(k) ~ k in the long 

wave-length limit, where k is the magnitude of the wave particle potential V(rt 



vector, while in simple liquids S(k) is quadratic with 
5(0) > 0. The vanishing of 5(0) indicates that jammed 
mono-disperse systems are hyperuniform [l2j]. 

To successfully avoid crystallization, a binary mixture 
of particles with two sizes (bi-disperse) is usually em- 
ployed. Although some mechanical and vibrational prop- 
erties of marginally jammed solids do not qualitatively 
depend on the ratio of the two particle sizes, we are 
short of systematic knowledge about how the particle- 
size ratio might affect jamming. For bi-disperse systems, 
it is also tempted to assume that either the large or the 
small particles can be taken as representative to describe 
the properties of the whole system. For instance, one 
may assume that S(k) measured for large particles and 
small particles are identical, so that the observation of 
S(k) ~ k in jammed mono-disperse systems is universal 
and all jammed systems are thus hyperuniform. 

In this paper, we study whether and how the particle- 
size ratio affects the properties of jammed binary mix- 
tures. By increasing the ratio R of the diameter of the 
two species from R = 1, we observe an increase of the 
critical volume fraction <f> c and a decrease of the bulk 
modulus accompanied with the increase of rattlers (par- 
ticles with no interacting neighbors). Most importantly, 
the inclusion of bi-dispersity qualitatively changes the 
static structure factor S(k) especially in the long wave- 
length limit. The linear behavior of S(k) in the long 
wave-length limit is thus only valid for jammed mono- 
disperse systems. 



II. METHOD 

The systems studied consist of N = 10, 000 frictionless 
spheres with the same mass m, contained in three dimen- 
sional cubic boxes with length L and periodic boundary 
conditions. Half of the particles have a diameter a, while 
the other half have a diameter Ra with R > 1. Particles 
interact via short range harmonic repulsion: the inter- 

: e (1 — Tij /<Jij) 2 /2 when the 
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separation between particles i and j, rtj is less than the 
sum of their radii cry = (o"j + <Jj)/2 1 and zero otherwise. 
This numerical model has been widely used in the studies 
of jamming, and is a typical model of granular materi- 
als, glasses, and foams. We vary the diameter ratio R 
from 1 to 4 to carry out a systematic study of its effects 
on the properties of the jammed systems. We note that 
most of the reported studies have focused on two specific 
values of R: 1.0 and 1.4. To generate jammed packings 
of frictionless spheres, we quickly quench random config- 
urations at high temperatures and fixed volume fraction 
to the local energy minima using L-BFGS minimization 
routine We tune the volume fraction until a pack- 
ing with desired pressure is obtained. In order to reduce 
quench rate dependence, we start at a volume fraction 
which is pretty close to the destination, so that tuning 
the volume fraction does not lead to significant particle 
rearrangements. We use a, m, and e as the unit of length, 
mass, and energy. 

For each R and desired pressure, we generate 
100 jammed configurations from independent initial 
conditions. The pressure is measured from P = 
-^ttT! ' ^nT^-, where the sum is over all interacting 
particle pairs. We obtain the average value of the vol- 
ume fraction cf> and coordination number z over these 
100 configurations at fixed pressure. The static structure 

factor is calculated from S(k) = 




Ej-exp(i*'^) 



FIG. 1: (a) Rescaled pressure P/Po and (e) coordination num- 
ber beyond isostaticity z — z c versus compression above Point 
J cj> — 4> c for marginally jammed solids with the particle size 
ratio R = 1.0 (black circles), 1.2 (red squares), 1.4 (blue di- 
amonds), 1.6 (green up triangles), 1.8 (orange left triangles), 
2.0 (magenta down triangles), 2.5 (purple pluses), 3.0 (ma- 
roon crosses), and 4.0 (cyan stars). The black solid lines in 
(a) and (e) show P/P = (t^-^) 1 ' 2 and z-z c = 11(^-<^ C ) ' 6 , 
respectively, (b), (c), and (d) show the critical volume frac- 
tion at Point J <j> c , concentration of rattlers at Point J cq c , 
and Pq as a function of R, respectively. The red solid line 
in (d) show Pq oc 1 — co c - Po is in the units of e/o" 3 , where 



and averaged over configurations, where the sum is over a = (1 + R )a /2 is the average volume of particles 
all particles, and r 3 is the Cartesian coordinate of par- 
ticle j. Due to the periodic boundary conditions, the 
wave vector must be chosen as k = (n x x + n y y + n z z) 

k 



with the magnitude k = y n 2 , + + n 2 z , where 

n x ,n yi n z = 0, ±1,±2, .... Because of the particle 
bi-dispersity, we also separately measure the static 
structure factor 
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large-large particles S , Li(fc) = 
small-small particles Sss(k) = 
and large-small particles S^sik) — 



^-Re 



-ik 



where E; an d 



E; ex P (ifc ' n) E s ex P (" 
E s are over all large and all small particles, respectively. 
It is obvious that S(k) = [S LL (k) + S SS {k) + S LS {k)} /2. 



III. RESULTS AND DISCUSSION 

Figure [TJa) shows how the pressure P of marginally 
jammed solids with different R varies with the volume 
fraction <f> on the approach of unjamming transition. The 
critical volume fraction <p c is a fitting parameter to satisfy 
P = Pq (4> — 4> c ) . For systems with harmonic repulsion, 
it is widely accepted that a = 1 However, we have 
noticed in previous studies that a deviates from 1 when 
<t> — 4> c is not tiny. The pressure grows nonlinearly at 
high compressions, which leads to a w 1.2 in the range 
of volume fractions studied, as denoted by the solid line 



in Fig. [TJa). 

The value of R does affect the jamming transition and 
material properties of marginally jammed solids. When 
R increases from 1, <p c increases and reaches a maximum 
at R e (2.0,2.5), as shown in Fig.rjTb). When R -> oo, 
the total volume of small particles is zero, so the sys- 
tem will effectively be mono-disperse and consequently 
<^>c(oo) ~ ^c(l)- It is thus unsurprising that 4> C (R) is 
peaked at an intermediate R. By varying R, we obtain 
about 3.5% increase of <f> c from mono-disperse systems. 
Since Point J controls properties of jammed solids [HI 
and probably dynamics of supercooled liquids approach- 
ing the glass transition [16j |, it is worthwhile to investi- 
gate further how R affects the dynamics of supercooled 
colloidal suspensions, e.g. glass transition and equation 
of state of hard sphere colloids. 

In marginally jammed solids, there are inevitably a 
fraction of particles that are rattlers. The concentra- 
tion of rattlers cq increases with decreasing the volume 
fraction approaching <j> e . For mono-disperse systems, it 
has been claimed that the rattlers are necessary to keep 
jammed packings of hard spheres saturated [12j. When 
R is away from 1, we would expect that the concentration 
of rattles rises, because it is easier for small particles to 
stay in the vacancies formed by large network particles. 
We measure cq(4>) for each R studied and extrapolate 
the concentration of rattlers at Point J, cq c , as shown in 
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FIG. 3: Pair correlation function g(r) for jammed systems 
with R = 1 (bottom, black) and R = 1.4 at P = 1CT 4 . For 
R — 1.4, g(r) is measured for both large particles (top, red) 
and small particles (middle, blue), r is scaled with the particle 
diameter. g(r) for R = 1.4 systems are shifted up in order to 
distinguish the curves. 



FIG. 2: Concentration of particles with Z interacting neigh- 
bors of (a) all particles measured for R — 1 (black circles) 
and R = 1.4 (red squares) at P = 10 -4 , and of (b) large par- 
ticles (black diamonds) and small particles (red triangles) for 
R = 1.4 systems shown in (a). The lines are guide to eye. 

Fig. [TJc). cq c increases with P and finally saturates at 
co c = 0.5 when P is large, indicating that almost all the 
small particles are rattlers. When R ^> 1, large parti- 
cles tend to form jammed structure, while leaving small 
particles to stay anywhere in vacancies. It is thus triv- 
ial to study jamming at large R. From Fig. [TJc), we 
can also tell that the slope of co c (P) starts to decrease 
with increasing R approximately when 4> C {R) reaches the 
maximum. Therefore, the decrease of 4> c at large R in- 
dicates the trivial rattler's influence: the jammed solids 
are more mono-disperse like and less saturated. In the 
following, we will mostly concentrate on systems with R 
before <j) c (R) reaches the maximum. Even though the ef- 
fects of rattlers on jamming become trivial at large P, it 
is still interesting to know how the small particles affect 
dynamics and rheology of amorphous systems approach- 
ing jamming, since large particles are easier than small 
particles to jam. 

The rise of cq c with increasing R leads to a decrease 
of the pressure with increasing P, since fewer particles 
are involved in interactions. It is thus natural to imag- 
ine that the bulk modulus of marginally jammed solids, 
which is quantified by P , decreases with increasing R. 
Figure D^d) shows that Pq(R) decreases with increasing 
R while L is kept fixed. We express Pq in units of e/a 3 , 
where a 3 is the average particle volume, which is equiv- 
alent to fixing L or particle number density N/L 3 when 
varying R at fixed <f> and N . In the R — > oo limit, Po 
saturates at P (l)/2, because the system is effectively 
mono-disperse with only N/2 network particles while the 
volume of the system is kept fixed. Thus, Po(R) is solely 
determined by the concentration of rattlers. As shown 
in Fig.[ljd), Pq(P) is indeed proportional to 1 — cq c . We 



note that one can vary R while keeping the diameter of 
the small particles a fixed. In this case, L increases as 
R is increased and this increase in L leads to a decrease 
of P so that Po would decrease much faster with R and 
eventually becomes zero in the R — > oo limit. Special 
care is thus required to specify what length scale, L or 
a or even Per, is being kept fixed when one studies how 
the material properties such as bulk modulus and shear 
modulus vary with R. 

In Fig. [He), we show the coordination number z as 
a function of (f> — <f> c for the same systems in Fig. [TJa). 
Similar to what has been observed, z = z c + zo(<fi — <fi c ) b , 
but b w 0.6 which is slightly different from 0.5 as previ- 
ously reported for systems with harmonic repulsion [l|. 
z{(f)—4> c ) does not show any dependence on P. Therefore, 
the geometric properties represented by z is universal for 
all P. It has been proposed that z plays a key role in de- 
termining vibrational properties of jammed solids. The 
collapse of z(<f) — <f> c ) over a wide range of P enhances 
our expectation that the coordination number is a rele- 
vant order parameter to describe the jamming transition. 
From the scalings of P and z with <fi — <j> c , we can eas- 
ily tell that z oc P 1 / 2 . Therefore, the pressure should 
be a better tunable parameter than the volume fraction 
to characterize properties of marginally jammed solids at 
fixed P. 

The coordination number discussed above is a mean 
field value. In marginally jammed solids, the number 
of interacting neighbors per particle Z fluctuates over 
particles. Figure d^a) shows the concentration of parti- 
cles with Z interacting neighbors, cz for jammed systems 
with P = 1 and 1.4. Z spans over a significant amount 
of values. For a d— dimensional system, a sphere needs 
at least d + 1 constraints to be stable. In Fig. (Hb), we 
separately plot cz for small and large particles. Appar- 
ently, large particles tend to have more constraints than 
small ones, because they have larger surface area to ac- 
commodate more small particles. For the systems shown 
in Fig. (2fb) , small particles have an average coordination 
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number of 5.15, which is much lower than the isostatic 
value z c . Therefore, the less constrained small particles 
should be easier than the over constrained large particles 
to rearrange during plastic events when the system is 
under strain. From Fig.[2]Jb) it is obvious that small par- 
ticles contribute to the majority of rattlers, as discussed 
above. 

The large deviation of Cz between large and small 
particles implies possible structural difference between 
them. Figure [3] compares the pair correlation function 
g(r) for the same systems in Fig. [5] In mono-disperse 
systems, the second peak of g{r) splits into two subpeaks 
at r = \/3 and 2 whose right hand side drops discontinu- 
ously near Point J [l7|,ll8[, indicating the glass formation. 
In bi-disperse systems, g(r) for large particles looks sim- 
ilar to mono-disperse systems, except for several more 
subpeaks at r < V3 due to the mixing of particles with 
different sizes. However, g(r) for small particles is dif- 
ferent. The subpeaks at r < 2 shift to the right, so it is 
less probable for four small particles to form the cluster 
corresponding to the subpeak at r = 

The apparent variation seen in g(r) reflects that the 
counterpart of g(r) in the Fourier space, the static 
structure factor S(k), may vary with species and R 
as well. In Fig. [4] we plot the structure factor mea- 
sured for large-large particles Sll^), small-small parti- 
cles Sss(k), large-small particles S^sik), and all parti- 
cles S(k). Contrary to what one might have expected, 
Shhik) and Sss(k) are distinct, especially in the long 
wave-length limit. There are also strong dependences of 
R on all the structure factors. When the pressure is suffi- 
ciently low, the structure factors do not have observable 
pressure dependence. Therefore, we compare the struc- 
ture factors for different R at the same pressure P = 10~ 4 
without losing generality. 

In the range of R that we focus on (R £ (1, 2]), SLh{k) 
is liquid like in the long wave-length limit: Sll (fc) ss 
S L L(0) + {(k) 2 , as shown in Fig. g^a). Both Sll(0) and C 
decreases with increasing R. In thermal equilibrium sys- 
tems, the structure factor at k = is proportional to the 
isothermal compressibility. Although marginally jammed 
solids are out of equilibrium, 5z,l(0) still provides us with 
the information about how the large particles are spa- 
tially distributed. When R increases, large particles get 
closer, which leads to the decrease of the compressibility. 
When R is so large that almost all the small particles are 
rattlers, Sll(O) ~ is recovered and the jammed system 
is trivially equivalent to mono-disperse. 

In the long wave-length limit, Sss(k) presented in 
Fig. 0Jb) has completely different k dependence from 
Shhik). Sss(k) decreases with increasing k and 5*ss(0) 
increases with R. Sss(k) at low fc's can be approxi- 
mately fitted with the Ornstein-Zernike form, Sss(k) = 

Sss(0)/ 1 + (£fc) 2 ; where £ increases with R. When R 
increases, viewed from small particles, the system is more 
dilute due to the growing amount of rattlers, which leads 
to an increasing compressibility and Sss(®)- The length 
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FIG. 4: Structure factor for (a) large-large particles SLh{k), 
(b) small-small particles Sss(k), (c) large-small particles 
SLs{k), and (d) all particles S(k) measured at P — 10 -4 and 
at R = 1.02 (black solid), 1.05 (red solid), 1.1 (blue solid), 
1.2 (green solid), 1.4 (orange solid), 1.6 (magenta solid), 1.8 
(purple solid), and 2.0 (cyan solid). The maroon dashed line 
in (d) is for R — 1. The dot-dashed line in (d) has a slope of 
1. 



£ extracted from the Ornstein-Zernike approximation of 
Sss(k) implies that small particles tend to form clusters 
with an average length scale £ that slightly increases with 
R. Such a length scale would saturate to the size of the 
vacancies formed by large particles in the large R limit. 

From comparisons between Fig.[5fa) and (b), when R 
is close to 1, SLLik) and Sss(k) which are self-part den- 
sity correlations are almost identical. In the long wave- 
length limit, they are nearly constant in k. When R *~ 1, 
Sls(0) is negative and cancels Sll(0) + 5ss(0) to give 
S(k) ss 0. As shown in Fig. [U^c), Sls(Q) increases with 
R but remains negative in the range of R studied. 

The strong R dependence of Sll^), Sss(k), and 
SLs(k) indicates that what we have learned from mono- 
disperse jammed systems, e.g. S(k) ~ k, is not gen- 
eral and thus not special features of marginally jammed 
solids. To emphasize this point, we plot the structure fac- 
tor of all the particles S(k) in Fig. 0Jd). S(0) grows up 
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with increasing R, which is consistent with the change 
of the bulk modulus with R observed in Fig. [TJ The 
long wave-length behavior of S(k) also strongly depends 
on R, which undergoes the transition from an increasing 
function to a decreasing function of k. It is obvious that 
S(k) <~ k is correct only in the R—>1 limit. 

IV. CONCLUSIONS 

In conclusion, by tuning the particle-size ratio R, we 
are able to study the possible effects on the properties of 
marginally jammed solids. The critical volume fraction 
of the jamming transition is a nontrivial functions of R 
when R is not large. Most strikingly the density fluc- 
tuations strongly depend on R and the particle species. 
Only jammed mono-disperse systems (R = 1) are hyper- 
uniform with S(k) ~ k. 

Our observations of the material property and struc- 
tural change of the jammed solids with R have some in- 



teresting implications, especially in the studies of the dy- 
namics and rheology of amorphous systems on the ap- 
proach to jamming. Since the small and large particles 
have distinct density fluctuation and such a difference 
increases with R, it would be interesting to see the co- 
existence of jammed large particles and liquid- like small 
particle clusters under thermal perturbations. As an in- 
dependent variable, R may be one more parameter to 
control the scaling collapse of the relaxation time of glass- 
forming liquids [l6|. It is also important to know whether 
R has an effect on the equation of state, glass transition, 
and glass fragility of glass-forming colloidal suspensions 
and their response to the shear stress. 
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